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2 ■ Abstract 

, Higher-dimensional orthogonal packing problems have a wide range of practical 

I applications, including packing, cutting, and scheduling. Previous efforts for exact al- 

O ■ gorithms have been unable to avoid structural problems that appear for instances in 

^ . two- or higher-dimensional space. We present a new approach for modeling packings, 

I using a graph-theoretical characterization of feasible packings. Our characterization 

■ allows it to deal with classes of packings that share a certain combinatorial structure, 



instead of having to consider one packing at a time. In addition, we can make use of 
elegant algorithmic properties of certain classes of graphs. This allows our characteri- 
zation to be the basis for a successful branch-and-bound framework. 

This is the first in a series of papers describing new approaches to higher-dimensional 
packing. 
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1 Introduction 



The problem of cutting a rectangle into smaller rectangular pieces of given sizes is known 
as the two-dimensional cutting stock problem. It arises in many industries, where steel, 
glass, wood, or textile materials are cut, but it also occurs in less obvious contexts, such as 
machine scheduling or optimizing the layout of advertisements in newspapers. The three- 
dimensional problem is important for practical applications such as container loading or 
scheduling with partitionable resources. It can be thought of as packing boxes into a con- 
tainer. When arranging these boxes in a container C, we have to preserve the orientations of 
the boxes; this constraint usually arises from considerations for stability in packings ("this 
side up"), from asymmetric texture of material in cutting-stock problems, or from different 
types of coordinates in scheduling problems. In the context of technical computer science. 



Teich, Fekete, and Schepers (2001) consider an application of three-dimensional packing to 
dynamic reconfiguration of hardware; one of the axes corresponds to time, the two others to 
different spatial dimensions. 

We refer to the generalized problem in c? > 2 dimensions as the d- dimensional orthogonal 
knapsack problem (OKP-d). Being a generalization of the bin packing problem (BPP), the 
OKP-fi is A/'P-complete in the strict sense - see Garey and Johnson (1979)j The vast ma- 
jority of work done in this field refers to a restricted problem, where only so-called guillotine 
patterns are permitted. This constraint arises from certain industrial cutting applications: 
guillotine patterns are those packings that can be generated by applying a sequence of edge- 
to-edge cuts; see [Christofides and Whitlock (1977)| and |Wang (1983)1 The recursive struc- 
ture of these patterns makes this variant much easier to solve than general or non-guillotine 
problems. 

Relatively few authors have dealt with the exact solution of non-guillotine problems. All 
of them focus on the problem in two dimensions. Biro and Boros (1984) [ gave a characteriza- 
tion of non-guillotine patterns using network flows but derived no algorithm. Dowsland (1987)] 



proposed an exact algorithm for the case that all boxes have equal size. [Arenales and Morabito (1995)| 
extended an approach for guillotine problems to cover a certain type of non-guillotine pat- 
terns. So far, only two exact algorithms have been proposed and tested for the general 
case of knapsack problems. Beasley (1985)| and Hadjiconstantinou and Christofides (1995)| 
have given different 0-1 integer programming formulations for this problem. Even for small 
problem instances, they have to consider very large 0-1 programs, because the number of 
variables depends on the size of the container that is to be packed. The largest instance that 
was solved in either article has 9 out of 22 boxes packed into a 30 x 30 container. (See the 
"OR-Library" by Beasley (1990)| for a number of test instances.) After an initial reduction 
phase, Beasley (1985) gets a 0-1 program with over 8000 variables and more than 800 con- 
straints; the program by Hadjiconstantinou and Christofides still contains more than 1400 
0-1 variables and over 5000 constraints. From Lagrangian relaxations, they derive upper 
bounds for a branch-and-bound algorithm that are improved using subgradient optimiza- 
tion. The process of traversing the search tree corresponds to the iterative generation of an 
optimal packing. 

We should note that more recently, two papers on the related problem of two- and three- 
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dimensional bin packing have been presented: Martello and Vigo (1998)] consider the two- 
dimensional case, while Martello, Pisinger, and Vigo (2000) deal with three-dimensional bin 
packing. We discuss aspects of those papers in ( Fekete and Schepers 1997c|, when consid- 
ering bounds for higher-dimensional packing problems. Even more recently, [Padberg (2000)] 
gave a mixed integer programming formulation for three-dimensional packing problems, sim- 
ilar to the one anticipated by the second author in his thesis ( Schepers 1997| ). Padberg 
expressed the hope that using a number of techniques from branch- and- cut will be useful; 
however, he did not provide any practical results to support this hope. 

In this paper we describe a different approach to characterizing feasible packings and 
constructing optimal solutions. We use a graph-theoretic characterization of the relative 
position of the boxes in a feasible packing. This allows a much more efficient way to construct 
an optimal solution for a problem instance: combined with good heuristics for dismissing 
infeasible subsets of boxes, we develop a two-level tree search. This exact algorithm has been 
implemented; it outperforms previous methods by a wide margin. 

The rest of this paper is organized as follows: after a formal problem description in 
Section |21 we introduce the novel concept of packing classes in Section El and show that it 
suffices to deal with the existence of feasible packings. In Section El we describe algorithmi- 
cally useful properties of packing classes, and highlight algorithmic difficulties in Section 
In dFekete and Schepers 1997dl ), we give a detailed description of how packing classes can be 
combined with other algorithmic ideas for obtaining powerful exact algorithms for orthogonal 
packing problems. 



2 Preliminaries 

Why are two- and higher- dimensional packing problems harder to solve than their one- 
dimensional counterparts? 

In order to decide whether another item can be packed into a partially filled container, 
in the one-dimensional case, we only need to know the total size of the items that have 
already been placed. When dealing with a two- or higher-dimensional problem, we also 
need to know the arrangement of boxes, i.e., their packing. Thus, the performance of an 
algorithm for higher-dimensional packing hinges on the use of an appropriate representation 
of packings. 

A main difficulty is highlighted by the following easy observation: In general, the feasible 
space for packing further objects into a container is not convex. This means that we cannot 
simply formulate a higher-dimensional orthogonal packing problem as a straightforward inte- 
ger linear program. Nevertheless, Beasley (1985)| and Hadjiconstantinou and Christofides (1995")] 



have used such a formulation, where the free space is controlled with the help of a large num- 
ber of 0,1- variables, describing a partition of the feasible space into convex pieces. They only 
consider the two-dimensional case. These formulations lead to very large 0,1-programs, even 
for packing problems of moderate size. It seems hopeless to try this approach for three- 
dimensional instances of relevant size. 

We present an alternative approach to these 0,1-formulations. Our new way of modeling 
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packings is based on a graph-theoretic characterization of packings. It can be used for any 
number of dimensions and leads to the central term of packing class, which describes a whole 
family of combinatorially equivalent feasible packings. The implementation of the resulting 
algorithms is described in detail in ( |Fekete and Schepers 19975| ). 
We start by introducing some basic notation and definitions. 



2.1 Problem Description and Notation 

In the following, we consider a set of items that need to be packed into a container. We 
concentrate on (i-dimensional boxes, and the items are also called boxes. 

The input data for a (i-dimensional orthogonal packing problem is a finite set of boxes V, 
and a (vector-valued) size function w : V ^ Qq that describes the size of each box in any 
dimension Xi, . . . ,a;„. For the orthogonal knapsack problem (OKP), we also have a value 
function v : V ^ that describes the objective function value for each box. 

The size of the container is given by a vector W G Q"*" . Without loss of generality, we 
assume that each individual box fits into the container, i. e., w{b) < W holds for each box. 

For the volumes of boxes and container we use the following notation. If b & V, and w is 
a size function defined on V, then 



voko{b) := Y\_^iib) 



i=l 



denotes the volume of box b with respect to w. Similarly, the volume of the container is 
denoted by 

d 

volw ■■= Yl^' 

i=l 

If 5 is a finite set and / a real-valued function on S, then we use the abbreviation 

f{S) :=5^/(x). 



xes 



Finally, we use the notation {v,w} for an edge in an undirected graph; we write {v,w) 
when referring to a directed edge. For a given set of vertices S ^ V in a graph G = {V, E), 
an induced subgraph is given by G[S] = {S, -^[5*]), with E[S] being the set of edges in G that 
connect vertices in S. For a vertex t> G ^ in a directed graph G = (V, E), S~{v) denotes the 
indegree of v. More graph-theoretic terms will be introduced when needed. In all cases, see 
the book by Golumbic (1980)] for more background. 



2.2 Orthogonal Packings 

We consider arrangements of boxes that satisfy the following constraints: 
1. Orthogonality: Each face of a box is parallel to a face of the container. 
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2. Closedness: No box may exceed the boundaries of the container. 



3. Disjointness: No two boxes may overlap. 

4. Fixed Orientations: The boxes must not be rotated. 

In the following, we imply these conditions when "packing boxes into a container" , "consid- 
ering a set of boxes that fits into a container" , and speak of packings. 

In order to formalize the notion of a packing, we introduce a Cartesian coordinate system 
with axes parallel to the edges of the container. The origin is one corner of the container, 
which is fully contained in the positive orthant. Thus, a packing is a function that maps 
each box to the coordinates of the one of its corners that is closest to the origin. 

Definition 1 (Packing) 

Given a set of boxes V, a size function w, and a container size W. 
A function p : V ^ (Qq is called a packing of (V, w, W), iff 

ybeV: p{b)+w{b)<W (1) 
yb,ceV,b^c,3te{l,...,d}: If{b) n if (a) = 0. (2) 

Here, If{b) denotes the interval [pi{b),Pi{b) + Wi{b)), for a box b ^ V, and a direction 
ie{l,...,d}. 

Condition ((T]) implies closedness, while (j2I) implies disjointness of the boxes. Because 
components of the size function are fixed, we pack with fixed orientations. 

Remark 1 The condition "p is a packing for (V, w, W) " can be expressed by linear con- 
straints. For this purpose, condition ^ has to be replaced by 

d 

yb,ceV,b^c\J p,{b)+w,{b) <pi{c) V p,{c)+w,{c) <pi{b). 
1=1 



In Nemhauser and Wolsey (1988)[ p. 12f., it is described how the disjunction can be trans- 



formed into 2d + 1 linear inequalities by using 2d additional 0,1-variables. The resulting 
mixed integer programs seem to be beyond practical size. Our own experience with a mixed 
integer approach has been far from promising, see [Schepers (1997)1 it remains to be seen 



whether the more recent revival by Padberg (2000) will lead to any computational results. 



2.3 Objective Functions 

Depending on the objective function, we distinguish three types of orthogonal packing prob- 
lems: 

Definition 2 (Optimization Problems) 
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The Strip Packing Problem (SPP) asks for the minimal height Wd of a container 
that can hold all boxes, where the size in the other d — 1 dimensions Wi, . . . , Wd-i are 
fixed. 

For the Orthogonal Bin Packing Problem (OBPP), we have to determine the 
minimal number of identical containers that are required to pack all the boxes. 

In the Orthogonal Knapsack Problem (OKP), each box has an objective value. A 
container has to be packed such that the total value of the packed boxes is maximized. 



This classification follows the scheme introduced by Wottawa (1996)| (He introduces a 



fourth class called Pallet Loading, where all boxes have identical size and value 1; this is a 
special case of the knapsack problem.) To clarify the dimension of a problem, we may speak 
of OKP-2, OKP-3, OKP-rf, etc. 

Problems SPP and OBPP are closely related. For d G IN, an OBPP-(i instance can be 
transformed into a special type of SPP-((i + 1) instance, by assigning the same x^+i-size to 
all boxes. This is of some importance for deriving relaxations and lower bounds. 

For all orthogonal packing problems we have to satisfy the constraint that a given set of 
boxes fits into the container. This subproblem is of crucial importance for our approach. 

Definition 3 (Orthogonal Packing Problem (OPP-c/)) 

For a set of boxes V , a size function w, and the container size W , is there a packing for 

iy, w, w) ? 

Clearly, the GPP is the decision problem for the above optimization problems. 



3 Problem formulation and mathematical approach 
3.1 Modeling the Problem 

In the following, we describe a new way of modeling feasible packings. The basic idea is to 
use the combinatorial information induced by relative box positions. For this purpose, we 
consider projections along the different coordinate axes; overlap in these projections defines 
an interval graph for each coordinate. Using properties of interval graphs and additional 
conditions, we can tackle two fundamental problems of exact enumeration algorithms: 

1. How can we prove in reasonable time that a particular subset of boxes is infeasible for 
packing? 

2. How do we avoid treating equivalent cases more than once? 
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3.2 Packing classes 



Instead of dealing with single packings we handle classes of packings that share a certain 
combinatorial structure. This structure arises from the way different boxes in a packing 
can "see" each other orthogonal to one of the coordinate axes. (So-called box visibil- 
ity graphs have been considered as means for representing graphs, e.g., see the references 
in Fekete and Meijer (1999)] ) We will see that transitive orientations of certain classes of 
graphs correspond to possible packings if and only if specific additional conditions are sat- 
isfied. This allows us to make use of a number of powerful theorems (jGhouila-Houri 1962| 
IGilmore and Hoffmann 1964| IGolumbic 1980| Korte and Mohring 1989 1 to perform pruning 
of our branch-and-bound tree. 

For a d- dimensional packing, consider the projections of the boxes onto the d coordinate 
axes Xj. Each of these projections induces a graph Gi = {V,Ei): 

{6,c}Gi?i^/f(6)n/f(c)^0. 

In other words, two boxes are adjacent in Gi, if and only if their Xi projections overlap. (See 
Figure H for a two-dimensional example.) By definition, the Gi are interval graphs, thus they 
have algorithmically useful properties that are described in Section |31 
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parallel projections. 



A set of boxes C \/ is called Xi-feasible, if J2b€s'^i(^) — boxes in S can 

be lined up along the without exceeding the Xj-width of the container. Finally, it 

follows from the disjointness condition (j2)) that no two boxes may overlap in all coordinate 
projections. Thus, we get the following necessary conditions: 

Observation 1 For any feasible packing, the graphs Gi = {V,Ei),i = l,...,d have the 
following properties: 

PI : Each Gi := (V^, Ei) is an interval graph. 
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P2 : Each stable set S of Gi is Xi-feasible. 

d 

P3: fl^. = 0. 

i=l 

The main result of this paper is the converse: These three conditions are also sufficient. 

A family of E = {Ei, . . . , Ed) of edges sets for a vertex set V representing a set of boxes 
is called a packing class for {V,w), if and only if it satisfies the conditions PI, P2, P3. For 
any given Gi = (V, Ei), we denote by Gi the complement graph {V, Ei) for Gi. See FigureEl 
If Gi arises from a packing, any edge in Gi corresponds to two boxes with non-overlapping 
Xj-projection, hence we can think of Gi as the comparability graph of the boxes in direction 
Xi'. Two boxes are comparable (and hence there is an edge in G,), if their Xj-projections do 
not overlap, and one projection lies strictly to the left of the other projection. 

Next we consider directed graphs arising from some of these undirected graphs, indicating 
the partial order of the Xj-projections. A directed graph D = {V, A) is called a transitive 
orientation of a graph G = (V, E), iff there is a bijection O : E A with the property 

(a, b) e AA{b,c) e A^ (a, c) G A. 

Note that this implies that the undirected edge {a, c} must also be contained in E. For more 
background on transitive orientations, see the book by |Golumbic (1980)[ 

For any packing class E, we consider transitive orientations Fj of the comparability graph 
Gi, and call F = [Fi, . . . , Fa) a transitive orientation of the packing class E. In the following, 
we show that for any transitive orientation F of a packing class, there is a packing. For this 
purpose, define a mapping pf : V ^ IRg by 

pf{v) := max{pf(M) + Wi{u)\{u,v) G Fi} (3) 

for V eV, i e {1, . . .,d}. 

3.3 Sufficiency of Packing Classes 

Without loss of generality, it is sufficient to consider packings that are generalized "bottom- 
left" justified, such that any lower bounding coordinate of a box is at zero or at the upper 
bounding coordinate of another box (without requiring the two boxes to touch). More 
formally, we say that a packing is gapless, if 

"^i e {1, . . . ,d} \fv eV : pi{v)=0 V 3u e V : Pi{v) = Pi{u) + Wi{u) . 

In the following, we show that the above three conditions on the graphs forming a packing 
class are not only necessary for the existence of a packing, but they are also sufficient, even 
for the subclass of gapless packings. This implies that a packing class is not just a set of 
graphs, but can also be interpreted as a whole equivalence class of feasible packings. 

Lemma 1 Any packing class E represents a feasible packing. More precisely, for a transitive 
orientation F of E, the function p^ : V ^ JRq describes a gapless packing. 
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Proof: Consider a packing class E. As each Gi = {V, Ei) is an interval graph, the 
complement Gi is a comparability graph, hence it has a transitive orientation Fj. Let F = 
{Fi, . . . , Fd) be any such set of transitive orientations. As defined above, this induces a 
mapping p^. We place all boxes v at the positions given by p^{v). To show that is a 
packing, we have to show that 

1. All boxes are within the boundaries of the container. 

2. No two boxes overlap. 

To establish 1., we show that at any stage there must be a set of boxes that determines 
the overall Xj-width of the partial packing, so by P2, the partial packing must fit in all 
directions. For this purpose, let f G V, i G {1, . . . ,d}. By construction of p^, the acyclic 
digraph (V, Fj) contains a directed path {v^^\ . . . , v^'^^) with 6~{v^^^) = and v^'^^ = v, such 
that 

r-l 
k=0 

Because Fj is transitive, S := {f . . . , f ^''^} induces a clique in Gi, implying that S" is a 
stable set in Gi. Hence by condition P2, 

r 

pf i^) + '^ii^) = < Wi, implying 1. 

fe=0 

To see 2., consider u,v E V, u v. By P3, there is an z G {1, . . . , c?} with the (undirected) 
edge {u,v} ^ Ei, hence {u,v} G Fj. Because Fj is an orientation of Fj, either {u,v) G Fj or 
{v,u) G Fj. By construction of p^, we get 

pf (v) > Pi{u) + Wi{u) V p[{u) > p[{v) + Wi{v). 

In both cases, u and v can be separated by an Xj-orthogonal hyperplane. Hence, the boxes u 
and V are disjoint and p^ describes a packing. By construction of p^, the packing is gapless. 
□ 

From Observation ^ and Lemma ^ the main theorem of this paper follows: 

Theorem 1 A set of d- dimensional boxes V can be packed into a container, iff there is a 
packing class E for (V, w), i.e., a set of d graphs Gi, . . . ,Gd with the properties PI, P2, P3. 

The proof of Lemma Q shows constructively that any orientation of a packing class cor- 
responds to a packing; basically, an orientation (fi,f2) of an edge {^1,^2} G Fj means that 
f 1 can "see" V2 along the positive Xj-axis. To illustrate the structure, we give the following 
example: 
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Example 1 Consider the OPP-2 P = {V,w), defined by 

V := {bi,b2,b3,b4,h} 

wib,) := (4, 1) , wib2) := (5, 1) , wih,) := (1, 3) , 
wib,) := (2, 2) , wib,) := (1, 2) , W := (5, 5) . 

and the packing class E for P as shown in Figure^ (Note that box i is darker than box j 
whenever i < j .) 



(VM 

• I I • 




(V,A2) 

(a) Cb) (c) (d) 

Figure 2: Constructing a packing for a two-dimensional packing class E: (a) A pair of in- 
terval graphs, (b) The corresponding complement, i.e., comparability graphs, (c) Transitive 
orientations of the comparability graphs, (d) A packing built from the transitive orientations. 
It is feasible by properties P2 and P3. 



The complements of the component graphs Ei and E2 from the example in Figure\^ each 
have six transitive orientations. Hence there are 36 transitive orientations of E. Figure\^ 
shows the corresponding packings, constructed by virtue of 

Remark 2 Independent from our work, an approach with some similarity to our pack- 
ing classes can be found in an unpublished report by Jerrum (1987)\ In that report, two- 



dimensional packings are described as pairs of complements of partial orders. However, 
Jerrum does not observe that there is a purely combinatorial characterization that allows it 
to map graphs to packings, as described in our LemmaU[ The latter is the main result of 
this paper. 

We conclude this section by a graph-theoretic formulation of a useful geometric property. 
If a set of boxes S has a total sum of Xj-widths exceeding k times the Xj-width of the 
container, then any feasible packing must have an Xj-orthogonal cut that intersects at least 
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Figure 3: All feasible packings of the packing class E from Figure El (a) 



k + 1 boxes. Using the terminology of packing classes, we get an algorithmically useful 
formulation of this condition: 

Theorem 2 Let E be a packing class for {V, w), i & {1, . . . , d}, Gi := {V, Ei) and S V. 
Then Gi[S] contains a clique of size 



Proof: As an interval graph, Gi is perfect. Thus, the induced subgraph Gi[S] has 
maximum clique size co := u{Gi[S]) equal to the chromatic number x{Gi[S]). Hence, S can 
be partitioned into uj disjoint stable sets: S = U^^^S'^. All stable sets from Gi[S] are also 
stable in Gi. Then P2 implies for all k E {1, . . . that Wi{Sk) < 1. Altogether, we get 

(UJ \ LU 

[J Skj = ^Wd{Sk) < UJ. 
k=l J k=l 

Because u is integer, the claim follows. □ 
A more elaborate discussion on bounds for higher-dimensional packing and their appli- 
cations can be found in Fekete and Schepers (199fc)| 



4 Algorithmic Implications 

When trying to solve a packing problem to optimality, a crucial step is deciding the orthog- 
onal packing problem (OPP) whether a particular set of boxes can be packed in a feasible 
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way. In a practical context, one will typically use three different tools: 

(1) In order to see whether there is a fast positive answer we can try to find a packing 
by means of a heuristic. 

(2) A fast way to get a negative answer is described in Fekete and Schepers (199fc)j 
where we discuss lower bounds for packing problems. 

However, for hard instances, both of these approaches will fail to provide an answer, 
meaning that a third tool has to be applied: 

(3) If all else fails, we need to resort to an appropriate enumeration method to solve the 
OPP. 

In this section, we describe how packing classes can be used to provide (3) if both of the 
easy approaches fail. (Because the OPP is NP-hard in the strong sense ( Garey and Johnson (1979)] ), 
it is reasonable to use enumerative methods.) 

As we saw above, the existence of a packing is equivalent to the existence of a packing 
class. Furthermore, we have shown that a feasible packing can be constructed from a packing 
class in time that is linear in the number of edges. This allows us to search for a packing 
class, instead of a packing. As we describe in the following, the advantage of this approach 
lies not only in exploiting the symmetries arising from different transitive orientations of the 
same packing class, but also in the fact that the structural properties of packing classes give 
rise to very efficient rules for identifying irrelevant portions of the search tree. 

Some of the technical implementation details of the enumeration scheme are rather in- 
volved; they are described in ( |Fekete and Schepers 199'7d |. Here we give an overview over 



the underlying mathematical ideas. 



4.1 Building Packing Classes 

When trying to build a packing class Ei, . . . ,Ed for a set of boxes, we use a branch- and- 
bound approach. The algorithm branches by fixing any edge {b, c} G Ei or by excluding it 
via {b,c} ^ Ei. The set of edges that are "positively fixed" to be in some E^ is called E^, 
while the set of "negatively fixed" edges that are excluded are called E_. Fixing an edge to 
be in is called augmenting E+, while excluding it means augmenting E^. When referring 
to a particular E^, the respective sets of included and excluded edges are denoted by E^^i 
and E_^i. 

The branching has two possible objectives: 
(Y) Augment E^ until it is a packing class. 

(N) Prove that no augmentation of E^ to a packing class is possible, as it would have to 
use "excluded " edges from E_. 

In the first case, our tree search has been successful. In the second case, the search on the 
current subtree may be terminated, because the search space is empty. Otherwise, we have 
to continue branching until one of the two objectives is reached. 



12 



In the following we show how graph-theoretic properties can be exploited to help with 
these objectives. 



4.2 Excluded Induced Subgraphs 

Following the idea described in the previous paragraph, we need three components for our 
enumeration scheme: 

(1) A test "Is the current set of fixed edges a packing class?" 

(2) A sufficient criterion that has no feasible augmentation. 

(3) A construction method for feasible augmentations. 

All three of these components can be reduced to identifying or avoiding particular induced 
subgraphs. 

First of all, it is easy to determine all edges that are excluded by condition P3. By 
performing these augmentations of E^ immediately, we can guarantee that P3 is satisfied. 
Thus we can assume that P3 is satisfied, and will remain satisfied by further augmentations. 

P2 explicitly excludes certain induced subgraphs: i-infeasible stable sets, i. e., i-infeasible 
cliques in the complement of each component graph. 

In order to formulate PI in terms of excluded induced subgraphs, we use the following 
powerful Propositions ^ and |H - the reader is referred to the book by Golumbic (1980)] In 
this context, we use a couple of technical terms: 

Definition 4 (Chords, Cocomparability Graphs) 

For a cycle C := [6o, • • • , ^fc-i, bk = ^o] of length k, the edges bibj,i,j G {0, . . . , k — 1} with 
(l^— j| mod k) > 1 are called chords; the chords bib j,i,j G {0, . . . , k—1} with {\i—j\ mod k) = 
2 are called 2-choids of C . A cycle is (2-) chordless, iff it does not have any (2-) chords. 

A graph is called triangulated, if it does not contain any chordless cycle of length k > 4. 

A graph G = {V, E) is a cocomparability graph, if the complement graph G = {V, E) is 
a comparability graph. 

In a graph G = (y,E), {vi,V2,V3} C V is called an asteroidal triple, if any two of the 
vertices can be connected by a path that contains no vertex adjacent to the third vertex. 

Proposition 1 (Gilmore and Hoffman 1964) 

A cocomparability graph is an interval graph, iff it does not contain the chordless cycle C4 
of length 4 as an induced subgraph. 

Proposition 2 (Ghouila-Houri 1962, Gilmore and Hoffman 1964) 

A graph is a comparability graph, iff it does not contain a 2-chordless cycle of odd length. 

In the proof of Theorem |2l we will make use of yet another characterization of interval 
graphs. 

Proposition 3 (Lekkerkerker and Boland 1962) 

A graph is an interval graph, iff it is a triangulated graph that does not contain any asteroidal 
triple. 
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Thus, £'+ is a packing class, if for alH G {1, . . . ,d} the following holds (recall that P3 is 
assumed to be satisfied): 

1. (y, Ej^^i) does not contain a C4 as an induced subgraph. 

2. (V, does not contain an odd 2-chordless cycle. 

3. {y, E_^i) does not contain an i-infeasible clique. 
For illustration, we provide the following example: 

Example 2 See Figure^ Consider V = {f 1, f2, '^'3, "^4} and suppose that at some stage, we 
have 

= {{VU ^2}, {f3, V4}, {f4, Vi}} 

and 

E_^i = {{V2,V4}, {vi,V3}}. 

This leaves the decision for the edge {^2,^3}. Adding {^2,^3} to E+^i would create an induced 
C4 in {V,E^^i); by Theorem\^ this is not admissible, forcing {v2^V3} to be in 



7 




f ^3 



E . (excluded edges) 



. (included edges) 



V4 

Figure 4: Making use of Proposition [TJ The edge {^2,^3} must be added to E^^i to avoid 
an induced C4. 



In dFekete and Schepers 19975 1, we describe the technical details how these properties 
can be implemented for achieving a fast algorithm. The computational results described 
there show that the overall code outperforms previous approaches by a wide margin. 



5 A Complexity Result 

Tackling the OPP by building a packing class resolves some of the difficulties, but of course 
there is little hope of solving the problem in polynomial time by any method. We conclude 
this paper by pointing out the exact nature of the remaining difficulty. 

Because properties PI and P3 appear easy to deal with, it is natural to suspect that 
this difficulty lies in property P2. Thus, it seems reasonable to hope for an efficient method 
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that preserves properties PI and P3 when augmenting a "partial packing class": Given a 
d-tuple of edge sets that satisfies P3, can it be augmented such that properties PI and P3 
are satisfied? (Note that P2 would be inherited to all augmentations, because any stable set 
after an augmentation must have been a stable set before.) A polynomial algorithm for this 
task would reduce the computational difficulty of solving the OPP by means of enumeration. 
Formally, we have the following problem: 

Definition 5 (Disjoint rf-Interval Graph Completion Problem (DIGCP-rf) 

GIVEN: graphs Gi = (V, ...,Gd = {V, 8+,d) with flti = 0- 

QUESTION: For i G {1, . . . , c?}, is there a superset Ei for each S+^i, such that 

nti Ei = 0, and the graphs d = {V, Ei), . . . ,Gd = (V, Ea) are interval 

graphs ? 

In the dissertation of the second author ( [Schepers 1997| ) it is shown that the existence of 
a polynomial algorithm for this problem is highly unlikely: The OPP remains NP-hard, even 
if we have a (i-tuple of edge sets that satisfies P2 and P3. This indicates that it is indeed 
PI that makes life difficult. 

Theorem 3 For d>2, the Disjoint d-Interval Graph Gompletion Problem is NP-complete. 
Proof: Recognizing interval graphs can be done in linear time ( Korte and Mohring 1989 1 , 



so the problem is in NP. 

For showing that the problem is NP-hard, we give a reduction of 3SAT. In the following, 
all notation is consistent with ( |Garey and Johnson 1979 ). 



Constructing a DIGCP-2 instance from a 3SAT instance: 

Consider a Boolean expression $ in conjunctive normal form with variables Ui, . . . ,Un and 
clauses Ci, . . . , Cm- The following DIGCP-2 instance can be constructed in polynomial time: 

1. For each clause Cj define six mfenor vertices Xj^k,x'ji^, k E {1, . . . ,3} and six exterior 
vertices yj,k,yj^ky ^ ^ {l? • • • ^S, }. Connect these vertices according to Figure El The 
black edges form the set Cj^\ and the gray edges form the set Cj^\ Thus we have 

3 

Cj = {{Xj^iXj^2}, {Xj^lXj^s}, {Xj^2Xj,3}} U \^{-Cj,kyj,k-,Xj^ky'j^k-i^'j,kyj,k^-^'j,ky'j,k}-i 

k=l 

3 3 

k=l 1=1 
l=/=k 

2. For each variable Ui define the four vertices Zi, z'^, Zi, z[. If Ui is the /cth literal in clause 
Cj, then Zi, z[ are connected with the exterior vertices j/j^fc and y'j ^. of Cj. If Ui is the /cth 
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Figure 5: The edges used for representing clauses. 



literal of clause connected with Uj^k and y'^^^. Thus, we get the two edge 

sets 

:= [ziUj^k, Ziyj,k, Ziy'j ,,, zly'jj Ui is the fc-th literal of Cj} 

and 

Uf' = {ziyj,k, z'iVj^k, ZiVj^k^ Ky'j,k I is the /c-th literal of Cj] 
These edges are shown gray in Figure El The black edges form the set 

Now the DIGCP-2 instance I is given by 

m 3 n 

y ■= {xj,k,x'j^k,yj,k,yj,k} ^ [j{zi,z'i,Zi, 4} 

j=l k=l i=l 

m n 

j=l i=l 

j = l i=l i=l 

Constructing a solution for 3SAT from a solution of DIGCP-2: 

Let now E2) be a solution of /. For i G {1, . . . , n} let 

TRUE, Ziz[ ^ El 



Ui 



FALSE, Ziz[ e El 
16 



yj,k yj,k (Cj,k=u,) 

• o o o ■ ■ ■ o o 




Figure 6: The edges used for representing variables. 

In order to show that this assignment of variables satisfies $, we will repeatedly use 
Proposition 01 In particular, the proposition implies that an interval graph cannot contain 
an asteroidal triple or a C4 as an induced subgraph. 

If Ui = TRUE, hence Zi,z'i ^ Ei, then Ziz'^ G Ei must hold; Otherwise the edge set 
{zi,Zi, z[, z-} would induce a C4 in (V, 

Let Cj be a clause in $. In Figure El consider the subgraph induced in 2 by the interior 
vertices of Cj that has the edge set Cj . We note the following. 

('2) 

1. No vertex of the path (xj^i, x^- 3, Xj^'i) is adjacent to Xj^s by an edge in 

2. No vertex of the path (xj^i, x^- 2, Xj^s) is adjacent to Xj,2 by an edge in 

3. No vertex of the path {xj^2-, x'j^i, xj^i) is adjacent to xj^i by an edge in Cj 

Thus, the vertices xj^i, Xj,2, Xj^s form an asteroidal triple, unless one or several edges prevent 
this by introducing adjacencies. More precisely, E2 must contain at least one of the edges 

('2) 

Xj^ix'j i, Xj^2x'j 2, Xj,3x'j,3 in addition to Cj , because all other edges that could destroy the 
asteroidal triple are in 

Without loss of generality let Xj^ix^ ^ G E2 and hence Xj^ix^- ^ Ei. In order to prevent 
the four vertices {xj^i, x^j, y^j} from inducing a C4 in (y,Ei), Vj^iy'j i G Ei must hold. 
Using the same argument on £'2, VjAVj,! ^ E2 implies that the edge ZiZ^ (or Ziz'^, resp.) 
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belonging to the first literal Ui (or -Uj, resp.) must be in E2, hence not in Ei. Corresponding 
to our chosen truth assignment, we get Ui = TRUE (or Ui = FALSE, resp.) Therefore, clause 
Cj is satisfied for any j, and hence $. 

Constructing a solution for DIGCP-2 from a solution of 3SAT. 

Consider a truth assignment of the Ui that satisfies $. Without loss of generality, let the 
first literal in each clause be TRUE, and let i G {1, . . . ,n}. 
If Ui = TRUE, let 




Ui'-^^ U Ea with A 



{yj,k, y'j^kl "^i ^"th literal of Cj}. 



If Ui 



FALSE, let 




It is easy to see that the corresponding subgraphs are interval graphs: 
For j G {1, . . . , m}, let 






yj.i 
yj,i 



Figure 7: Representing {V, Cj ) by intervals. 
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Let now 



Xj^2 



X 



X 



j3 I 

- (2) 

Figure 8: Representing (V, Cj ) by intervals. 



i=l j=l 



m 

(1) , , I I ^ (1) 



n n m 



i=l i=l j=l 

{V, El) and {V, E2) are interval graphs, because all their connected components are. 
S+^i C El and ^ -E2 follows immediately from the construction. 

For j G {1, . . . , m}, let Ui be the first literal from clause cj. By assumption, this literal has 

the value TRUE, implying that the edge yj^iy'j i does not belong to f/j or t/j, resp. Therefore, 
Ci is disjoint from these edge sets. This excludes the only possible overlap of Ei and E2, 
implying EiH E2 = 0. Therefore, {Ei, E2) is a solution of the DIGCP-2 instance /. □ 



6 Conclusions 

In this paper, we have introduced a new way of characterizing higher- dimensional packings 
of a set of boxes into a container. Using graph-theoretic structures, this characterization has 
a number of nice algorithmic properties. We describe in (Fekete and Schepers 19975) how 



these properties can be used for implementing an exact algorithm that is able to solve two- 
and three-dimensional packing problems of relevant size. Our implementation makes exten- 
sive use of new classes of lower bounds for higher- dimensional packing problem; the math- 
ematical background is described in Fekete and Schepers (1997c)[ The overall algorithm 



appears to have convincing performance; this should make the mathematical structures and 
properties described in this paper interesting to a wide audience. 

Our graph-theoretic characterization has also turned out to be useful in practical ap- 
plications in which there are additional spatial constraints in some of the dimensions. See 
Teich, Fekete, and Schepers (2001)] for an application arising from dynamic reconfiguration 



of hardware. Making further use of order-theoretic structures, we have been able to extend 

this approach to handle order constraints. Details are described in Fekete, Kohler, and Teich (2001)( 
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